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ABSTRACT 


This thesis develops a methodology for updating the Navy’s Logistics Factors File, 
which has been neglected in recent years and requires updating. This study is limited to 
Repair Parts (Class IX of the Department of Defense Supply Class Codes) for the 
following four classes of ships: CVN-68 (Nimitz class) Aircraft Carriers, CG-47 
(Ticonderoga class) Guided Missile Cruisers, DD-963 (Spruance class) Destroyers, and 
FFG-7 (Oliver Hazard Perry class) Guided Missile Frigates. 

The current Logistics Factors File structure includes a single data entry in pounds 
per unit per day to describe the sustainment requirements of these units for all of the 
supply classes and their respective subclasses. For Repair Parts, these values are severely 
understated when compared to contemporary data. These “pounds per unit per day” 
random variables have heavily skewed distributions. These distributions can be fitted with 
mixtures of standard probability distributions, and it seems wise to recommend that 
associated variability information be included either directly in the Logistics Factors File, 


or in a readily available companion source. 
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EXECUTIVE SUMMARY 


This thesis develops a methodology for updating the Navy’s Logistics Factors 
File (LFF), which has been neglected in recent years and requires updating. The LFF 
contains logistics consumption planning factors for all units in the Navy by 
Class/Subclass of supply. It is used to predict unit sustainment requirements in 
support of Joint Operation Plans, crisis action planning, logistics feasibility 
estimation, and for wargaming analysis. This thesis is limited in scope to modeling 
only repair parts for the following four classes of ships: CVN-68 (Nimitz class) 
Aircraft Carriers, CG-47 (Ticonderoga class) Guided Missile Cruisers, DD-963 
(Spruance class) Destroyers, and FFG-7 Oliver Hazard Perry class) Guided Missile 
Frigates. Department of Defense guidance defines repair parts as Class IX, which is 
further subdivided into different subclasses. 

The current LFF structure includes a single data entry (in pounds per unit per 
day) to describe the sustainment requirements of these units for all of the Supply 
Classes and their respective Subclasses. For Repair Parts, these values are severely 
understated when compared to contemporary data. These “pounds per unit per day” 
random variables have severe positive asymmetries. These distributions can be fitted 
with mixtures of skewed probability distributions in order to achieve more reliable 
consumption figures. It is also wise to recommend associated variability information 
be included either directly in the LFF, or in a readily available companion source. 
The results obtained from using mixed distributions are exhibited alongside the 
current factors in the Table on the following page to illustrate the significant 


differences between the two: 








Factor Factor Factor Factor Factor Factor Factor Factor 


Current New Current New Current New Current New 
CVN-68 | CVN-68 CG-47 CG-47 | FFG-7 FFG-7 DD-963 DD-963 





Table Comparing Current Factors in Navy LFF To New Factors Derived re 
Recent Historical Data and Mixed Probability Distributions. All factors in 
pounds per unit per day. 


The methodology developed in this thesis fits mixtures of probability | 
distributions to historical data collected from the supply requisitions of Naval surface 
combatants over a three year period. Analysis performed on this data yields new 
factors which may be used to update the Navy’s LFF. The use of mixtures is 
recommended because of the substantial number of days where no requisitions occur, 
as well as the infrequent yet critical days where requisitions of ‘extreme’ weight 
occur. | 

There are several caveats that must be understood by potential users of these 
new factors.. These are explained in the development of the methodology and 
summarized in the conclusions. One of the main difficulties in the study was the 
discovery of missing weight data on many of the items that had been requisitioned 
during the data collection time frame. The other significant difficulty was the method 


of accounting for the extremely large requisition weights. Developing the methods 


XLi 








to correct these difficulties, as well as other, smaller ones, is a critical part of the 


methodology developed by this study. 
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I. INTRODUCTION 

This thesis fits mixtures of probability distributions to historical data taken from the 
supply requisitions of U.S. Navy surface combatants. The intent of this thesis is to devise a 
methodology for analyzing this data to update the Navy’s Logistics Factors File (LFF). The 
scope of this thesis is limited to modeling only repair parts for four common classes of 
warships. The use of mixtures is recommended because of the substantial numbers of days 
where no requisitions were submitted, as well as the infrequent yet critical days where 
requisitions of ‘extreme’ weight occurred. The results iad from using the mixed 
distributions are exhibited in Table 1 on the following page. However, there are many caveats 
that must be understood by potential users of the Table. These are explained in the 
development of the methodology and summarized in the conclusions. 

One of the main difficulties in the study was the discovery of missing weight on on 
many of the repair parts requisitioned during the data collection time frame. This problem, 
and how it was handled, is discussed in Section A of Chapter III. The other significant 
difficulty was the method of accounting for the outliers, or days with extremely large 
requisition weights. This issue is addressed in Section B of Chapter IV. 

Other aspects of the organization of this thesis are as follows: Chapter II contains 
general background information on the critical role of the LFF in Joint Operation Planning. 
More specific details about the Navy’s perspective with regards to the LFF and the 
Department of Defense (DoD) system of classifying Supply goods are contained in Chapter 
Tl. Data sitineliite manipulation, and the methods of analysis are in Chapter IV. These are 


followed by results of the study in Chapter V, and Conclusions in Chapter VI. Several 














appendices are included to document specific details of the data analysis phase. 













| Current 






Current | 
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Table 1. Goniparnon of Cunreat Factors i in Navy I LFF To New Factors Derived from 
Recent Historical Data and Mixed Probability Distributions. All factors in pounds per 
unit per day. 














Il. THE ROLE OF THE LFF ON JOINT OPERATION PLANNING . 


Logistics is defined as “the science of planning and carrying out the movement and 
maintenance of forces.” Military planners focus on logistics processes and products with the 
objective of achieving and sustaining operational readiness. In general, the mission of naval 
logistics is “to provide and sustain our operational readiness by getting the right support to 
the night place at the right time” [Ref 1]. 

The six functional areas of logistics are considered to be: supply, maintenance, 
transportation, erigineering, health services, and other services. Successful military planning 
involves careful integration of all six functional areas [Ref 2]. Of these six functional areas, 

Transportation has been a critical factor in strategy since fighting men carried 

equipment on their backs and lived off the countries where they were 

engaged. It has grown more important as the scope of hostilities has widened, 

and the burden of military equipment and supplies has increased [Ref 3]. 
Transportation resources are required to support mobilization, deployment, employment, 
sustainment, redeployment, and demobilization operations [Ref 4]. In short, no operation 
can succeed without transportation. 

A 1994 Joint Chiefs of Staff JCS) Deliberate Planning Conference concluded a need 
existed for the services to conduct a review of their LFFs. The LFF is one of several 
Standard Reference Files (SRFs) within the Joint Operation Planning and Execution System 
(JOPES), which supports planning and sustainment of forces. The LFF provides logistics 


consumption factors by class/subclass of supply for each unit in the Type Unit Characteristics 





(TUCHA) file, another JOPES SRF. (Table 2 on pages 11-12 lists the ten DoD supply 
classes and their respective subclasses.) These logistics consumption factors can be used in 
logistic feasibility estimation, wargaming analysis, or to produce a Time Phased Force and 
Deployment Data (TPFDD) database, which helps determine logistics sustainment 
requirements in the development of an Operation Plan (OPLAN ). The conclusion of the JCS 
conference was that existing factors were outdated and unreliable, invalidating logistics 
sustainability analyses [Ref 5]. Those familiar with DoD transportation systems and planning 
methods may choose to bypass the remainder of this chapter, which iat background on 


these issues and the importance of an accurate LFF. 


A. UNITED STATES TRANSPORTATION COMMAND 
The US Transportation Command (USTRANSCOM) is a single, unified cena 
whose mission is “to provide strategic air, land, and sea transportation to deploy, employ and 
sustain military forces to meet national security objectives across the range of military 
operations” [Ref 6]. Combatant commanders inform the Commander in Chief. U.S. 
Transportation Command (USCINCTRANS) of their movement requirements and required 
delivery dates (RDDs). USTRANSCOM then provides strategic lift through its three service 
component commands: Air Mobility Command (AMC), Military Sealift Command (MSC), 
and Military Traffic Management Command (MTMC). 
_ AMC, a responsibility of the Air Force, manages strategic airlift. U.S. war plans have 
traditionally required 5% of the highest priority cargo, mail and passengers to be airlifted into 


theater. Additionally, due to its speed and flexibility, airlift will be combined with pre- 














positioned assets to transport the initial U.S. forces deploying into a theater of conflict [Ref 
7|. The role of strategic airlift is limited by cost and capacity. For example, a typical strategic 
airlift mission from Operation Desert Shield/Storm consumed almost one million pounds of 
fuel and cost $280,000 [Ref 8]. The capacity limitations of strategic airlift are further 
constrained by the fact that only 20% of airlift capacity is contained in active duty forces. The 
other 80% is shared by reserve components (30%) and the commercial sector (50%), 
supported by the Civil Reserve Air Fleet (CRAF) [Ref 7]. 

MSC is managed by the Navy and provides strategic sealift, which is relatively 
inexpensive and has the highest capacity. Sealift has deaniioaally provided 95% of all 
strategic lift during force mobilizations. While the capabilities of airlift have improved, sealift 
still transported 90% of all materiel during Operation Desert Shield/ Storm [Ref 1]. However, 
despite vast upgrades to U.S. sealift capability during the 1980s, we still had to charter over 
70 foreign flagged vessels to support sealift requirements during Desert Shield. Much like 
airlift, only 20% of sealift capability is assigned to active duty forces. The other 80% is 
covered by the commercial sector [Ref 7]. 

The Army is responsible for MTMC, the third component of USTRANSCOM. 
MTMC manages ground transportation (truck and rail) within the Continental U.S. 
(CONUS), as well as most U.S. military ocean terminals [Ref 6]. In the event of a large scale 
mobilization of forces, MTMC assets would be involved in both the initial surge of troops and 


supplies as well as the follow-on sustainment of those forces. 








B. JOINT OPERATION PLANNING 

Joint operation planning involves the preparation and execution of plans for potential 
crises that involve military forces. The President of the United States directs the National 
Security Council (NSC) to establish national security procedures, including mobilization, to 
the DOD. From these directives, the Joint Planning and Execution Community (JPEC) plans 
for the mobilization, deployment, employment, sustainment, redeployment, and 
demobilization of joint forces [Ref 4]. 

Within the classification of joint operation planning are two classes of structured, 
formal processes: deliberate planning and crisis action planning. Crisis action planning is 
more reactive, usually carried out in response to a specific situation. Though crisis action 
planning develops in a very short time period, usually days or even hours, logistics 
requirements must still be met. Deliberate planning involves the entire JPEC community, and 
can take 18-24 months to fully develop [Ref 9]. The LFF is used in both methods of planning 


to validate transportation feasibility estimates and ensure forces can be adequately sustained. 


Joint 
Operation 
Planning 


Deliberate 
Planning 





Figure 1. Classes of Joint Operation Planning. 














USTRANSCOM conducts detailed transportation feasibility analysis during both 
deliberate and crisis action planning using the Joint Flow and Analysis System for 
Transportation (JFAST). JFAST is used during crisis action planning by providing 
information on transportation assets needed to move and resupply forces, as well as interpret 
shortfalls and assist in reallocating lift. Logistics planners only have to provide major units 
to be involved along with associated planning factors from the LFF. JFAST then uses its 

Notional Requirements Generator (NRG) to calculate combat support, combat service 
support, and non-unit cargo requirements. These requirements are then used to estimate the 
associated transportation requirements [Ref 10]. 

Fundamental to the deliberate planning process is the Joint Strategic Capabilities Plan 
(JSCP), which tasks the combatant Commanders-in-Chief (CINCs) wat preparing ae in 
support of national security objectives. There are four types of deliberate plans that can be 
prepared based on JSCP requirements. The type of plan developed is based on the level of 


detail required and purpose of the plan, as is shown in Figure 2 [Ref 11]. 


Deliberate 
Planning 


CONPLAN CONPLAN 
(Without Functional 
TPFDD) Plan 





Figure 2. Types of Deliberate Planning. 





The most complete and detailed joint operation plan is called an OPLAN. “OPLANS 
are normally prepared when the contingency has a compelling national interest and a specific 
threat, and the large scale of the contingency requires detailed prior planning for complex 
issues” [Ref 9]. Because OPLANs are so detailed and complex, the JOPES system which 
processes them requires many annexes and appendixes, including a TPFDD. A CONPLAN 
is a joint OPLAN in abbreviated concept format. Depending on the level of detail required 
by the JSCP, a CONPLAN may or may not contain a TPFDD [Ref 11]. The LFF plays a 
crucial role in the development of OPLANs and CONPLANSs (with TPFDDs), but is also used 
by JFAST to estimate lift requirements for CONPLANs (without TPFDDs) and Functional 
Plans. 

Ref 12, Joint Pub 5-03.1 “Joint Operation Planning and Execution System, Vol I’, 
defines a TPFDD as: | 

The JOPES data base portion of an operation plan; it contains time-phased 

force data, non-unit-related cargo and personnel data, and movement for the 

operation plan, including: 

@ In-place units. 

® Units to be deployed to support the operation plan with a priority 
indicating the desired sequence for their arrival at the port of 
débarkation. 

@ Routing of forces to be deployed. 


® Movement data associated with deploying forces. 


® Estimates of non-unit-related cargo and personnel movements to 
be conducted concurrently with the deployment of forces. 


e Estimate of transportation requirements by common-user lift 











resources, as well as those requirements that can be fulfilled by 
assigned or attached transportation resources. 

After planning movement requirements and the time-phased sequencing of forces into 
the theater, the TPFDD shows for the quantities of supplies that will be required to sustain 
those forces. These supplies are then also phased into the theater. The final phase in 
developing a TPFDD is validating the plan by simulating the force movements. In the event 
that no feasible strategic lift deployment can be derived, then shortfalls must be resolved by 
adjusting the movement requirements. The CINCs then determine whether such adjustments 
pose acceptable or unacceptable risks [Ref 9]. 

JOPES uses a modeling system called Logistics Sustainment Analysis and Feasibility 
Estimator (LOGSAFE) to provide feasibility estimates and generate resupply requirements. 
In doing so, LOGSAFE uses the LFF, TUCHA, Ports-of-Support File (POSF), Geographic 
Location File (GEOLOC), OPLAN TPFDDs and user-specified parameters [Ref 10]. Since 
ee om is such an important fundamental of logistics, and good logistics planning is 


essential to any military campaign, planners need an accurate LFF to plan for any contingency. 
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Ii. THE NAVY’S LOGISTICS FACTORS FILE 


The primary purpose of the LFF is to predict resupply planning factors for Navy and 
Marine Corps units. These factors must be accurate enough to ensure proper joint operation 
planning. The following list describes several reasons for maintaining an accurate LFF: 

A. Predicting sustainment requirements of non-unit related cargo for feasibility 
analysis of OPLANSs. 

B. Defining the lift footprint (.e., the impact of strategic lift capabilities) for 
determining amphibious lift requirements in support of future shipbuilding programs. 

C. Identify shortfalls in current OPLANs and reapportion supplies to correct them. 

D. Assess strategic lift requirements for apportioning limited sealift and airlift assets 
to commanders. 

E. Determine the required quantity and composition of prepositioned supplies to 
maximize slate 

F. Determine requirements for Host Nation Support (HNS) to meet identified 


shortfalls [Ref 10]. 


Table 2 lists the ten DoD classes of supply and their various subclasses [Ref 13]: 


CLASS SUBCLASSES 


I. Subsistence (food) A - Nonperishable dehydrated subsistence 
C - Combat rations requiring no organized dining facility 


R - Refrigerated subsistence 
S - Nonrefrigerated subsistence (less other subclasses) 
W - Water 





1] 











CLASS SUBCLASSES 


II. General support items: | A - Air 

Clothing, individual B - Ground support material 
equipment, hand tools, E - General supplies 
administrative, F - Clothing and textiles 
housekeeping supplies, etc. | G - Electronics 


M - Weapons 
T - Industrial supplies 


If. POL (Petroleum, Oils | A - Air (JP-5) 
& Lubricants): W - Ground (Surface) 
P - Packaged POL 


IV. Construction & barrier | A - Construction 
materials. B - Barrier materials 
V. Ammunition: A - Air 

Ammunition of all types W - Ground 


A - Non-packaged personal demand items 
M - Personal & Official letters and packaged mail 
P - Ration Supplement Sundry Packs (RSSP) 
















VI. Personal demand 
items: 
Non-military sales items. 






VII. Major end-items A - Air 
| B - Ground support material (generators, fire-fighting, etc) 
D - Admin & general purpose vehicles 


G - Electronics 
K - Tactical and special purpose vehicles 
L - Missiles 


M - Weapons 


VII. Medical material A - Medical material 
B - Blood and fluids 7 


IX. Repair parts A - Air 
B - Ground support material (generators, firefighting, etc) 
D - Administrative vehicles 
G - Electronics 
L - Missiles 
M - Weapons 
T - Industrial supplies 


NONE 


Table 2. Classes and Subclasses of Supply 











A. PROBLEMS WITH EXISTING LOGISTICS FACTORS FILE 

Despite its importance to operational planning, there are many shortfalls currently 
inherent to the Navy LFF. The following discussion of shortfalls was produced by the 
Potomac Research Corporation (PRC), Inc., in their second Interim Report for improving, 
validating and maintaining the LFF for the Navy [PRC, 1995]. 

1. Outdated Guidance 

The joint directive which provides the services specific guidance on providing inputs 
~ to the LFF is the Logistic Factors Report LOGFACREP), which was issued in 1985. The 
Chairman of the Joint Chiefs of Staff Manual (CJCSM) 3150.23 updates the LOGFACREP 
and recommends significant changes to the LFF. If these recommendations are approved, 
planning cycles will decrease from two years to one year, necessitating a simpler structure for 
the LFF, including the elimination of oe to “source” where non-unit related supplies cone 
from. 

2; Non-Standard Models for Generating Sustainment 

Since the late 1970s, records for TPFDDs have been developed by a system called the 
Movement Requirements Generator (MRG). In 1992, LOGSAFE was fielded as the 
replacement for the MRG and installed on the Worldwide Military Command and Control 
System (WWMCCS). While LOGSAFE is now the official model for generating TPFDD 
records, many planners on unified and service component command staffs still prefer to use 
MRG or a third model called Logistics Generator (LOGGEN), because they are easier to use 
or provide more details. This is permissible because of a loophole in the joint guidance which 


allows the use of other systems which provide similar data and analysis. 
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3; Lack of Information on the LFF 

The only published material that PRC, Inc., could find which provides any background 
information on the LFF was the old LOGFACREP. No information could be found 
describing the source of the intensity rates used by the LFF or the methodology employed 
(such as pounds per man per day for a given subclass of supply). Even the Joint Publications 
which describe the JOPES system (Joint Pubs 5-03.1 and 5-03.2) do not provide any of this 
background information. Hence, no insight could be determined as to why it is set i the way 
it 1S. 

4, No Routine Update Process for the Navy 

Original Navy logistic planning factors were produced from a 1979 study entitled “An 
Analysis of Navy Logistic Planning Factors”, dated 30 November 1979. Many of those 
factors still reside in the current Navy LFF. A partial update of the Navy LFF was sonducied 
in 1989 by the David Taylor Naval Research Center. Based on historical requisition data 
collected from Sepioned surface and aviation units, this appears to be the Navy’s most recent 
update. Despite the updated factors, this report was neither approved nor released, and none 
of the participants of that study could be contacted, so the exact methodology used by the 
study is unknown. 

5. Questionable Validity of Current Factors 

In addition to the current Navy factors being outdated and derived from unknown 
methods, their very low values question their accuracy. For example, the current factor for 
an FFG-7 Guided Missile Frigate, class IX(G) (Repair Parts, Electronics), moderate intensity, 


European conflict scenario is 14 Ibs per day. Historical data collected by the Ships Parts 
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Control Center (SPCC), Mechanicsburg, PA, reveals days where this type of unit 
requisitioned over 100,000 Lbs of this subclass. This magnitude of error is not limited to this 
class/subclass of supply, nor to this class of ship. 

6. Navy Uses Different Resupply System than Other Services 

It is important to note that the LFF is used to generate non-unit related cargo 
requirements. Non-unit related cargo is defined as “All equipment and supplies requiring 
transportation to an area of operations, other than those identified as the equipment or 
accompanying supplies of a specific unit” [Ref 12]. The TPFDDs generally address only the 
first 90 days of an OPLAN, during which most of the common user strategic lift provided by 
USTRANSCOM is dedicated to the Army (80%) and the Air Force (10%). These two 
services generally deploy with only three days of sustainment. Because of their need to be 
resupplied almost immediately after deploying, these services use a “push” system which 
provides resupply until a regular requisitioning system can be established. 

In contrast, Naval forces typically deploy with six months of sustainment carried by 
the combatant ships and the Combat Logistics Force (CLF) support ships. There are some 
notable exceptions to this policy, such as the Mobile Construction Battalions (MCBs) and 
Advanced Base Functional Components (ABFCs), but it applies to the Navy component of 
Amphibious Ready Groups (ARGs) and Carrier Battle Groups (CVBGs). This six months 
worth of sustainment is considered unit-related cargo because it accompanies the units. 
Hence, the Navy utilizes a “pull” resupply system, whereby it requisitions only what it needs 
as it operates. Naval forces still require sustainment, especially certain types of munitions, 


but this does not amount to a significant portion of common user strategic lift assets. 
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While this discussion may appear to lessen the importance of the Navy LFF with 
regards to joint operation planning, quite the opposite is true, especially for Class IX (repair 
parts). The Navy’s Class IX factors must be capable for use in the prediction of the 
sustainment requirements of repair parts that will be requisitioned for Battle Damage Repair 
(BDR). In such a case, such as a ship contacting a mine or being struck by a cruise missile, 
it is highly unlikely that it can be repaired using only assets carried by the Battle Group (BG). 
More likely, thousands of pounds of repair parts would have to be lifted into the theater to 
conduct the necessary repairs. Additionally, no ship or BG can be sietaihad to carry enough 
of every spare part to account for every conceivable cise casualty that a ship could 
suffer during a wartime deployment. The high priority that would be assigned to any repair 
part(s) needed to restore a ship to operational readiness would necessitate the part(s) being 
airlifted, and airlift assets will probably be in more demand than sealift. LFF sustainment ine 
are also used when planning logistics sustainment requirements for wargames. 

7. Navy Does Not Use DoD Class/Subclass Supply System 

Each subclass of supply listed in Table 2 has a consumption factor associated with it. 
It is difficult for the Navy to develop factors in this fashion because the Navy supply system 
does not assign items by class/subclass. Instead, the Navy has traditionally classified its 
sustainment supplies in the following broad categories: ordnance, fuel, dry cargo, wet cargo 
and refrigerated cargo. Navy supply infrastructure is decentralized to assorted System 
Commands (SYSCOMS), which assign personnel to manage certain assets. Because of this 
setup, the Navy cannot specifically manage supply by the DoD class/subclass system. The 


method that the Navy has devised to get around this problem is a critical part of the 
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methodology of this study, and it is discussed in Chapter IV. 

8. Inaccurate Supporting Files 

The LFF is used in conjunction with the TUCHA, which provides the Unit Type Code 
(UTC) for notional units. Each service is responsible for maintaining the data in their 
respective TUCHAs, but until recently the Navy has done a poor job of updating. PRC, Inc., 
found examples of UTCs with had consumption rates in the LFF that were not listed in the 
TUCHA, as well as UTCs in the TUCHA with no sustainment factors in the LFF. 
Additionally, the Navy TUCHA listed non-deployable units that would never require 
consumption rates, such as a Coast Guard District Headquarters. The Navy and Marine 
Corps even had different units in the TUCHA with the same UTC. Finally, as entire ship 
classes and aircraft squadrons are decommissioned, as well as new ship types are introduced 
(such as the DDG 51 and LHD 1 classes), a need exists to establish a system which 
periodically reviews and updates the TUCHA. 

Another area where incomplete data files hinder the Navy’s ability to maintain an 
accurate LFF is the lack of weight and cube (volume) information on items in the Navy and 
DoD supply inventories. The problems this causes will be discussed in detail in Chapter IV, 
but it should be apparent that it is impossible to develop consumption rates (in units of lbs per 
UTC per day) if the weights of the supply items are unknown. The Navy supply system listed 
469,361 items in its inventory (listed by National Item Identification Number, or NIIN). Only 
254,221 of these items (54.2%) had a weight assigned, while the other 215,140 had zeroes 
listed. The Defense Logistics Agency (DLA) did only slightly better, with actual weights 


assigned to 2,039,158 items out of 3,533,161 (57.7%). The remaining 1,494,003 items were 
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assigned weights of zero. 

9. No Consideration for Variance of Consumption Factors 

The factors listed in the LFF are point estimates of some central value, though it is 
unknown exactly what these numbers represent. They could be the mean or median values 
from some set of historical data, or even a best guess from someone with experience in this 
field. Regardless of how these numbers were derived, there is no consideration given to the 
very large variances that are involved. An example of the magnitude of this variance is cited 
in Subsection 5 of the present Chapter. 
B. THESIS OBJECTIVES 

The Navy is currently working to improve many of the problem areas discussed in 
Section A of this Chapter. Examples of these efforts include: 

@ Reviewing and providing comments on the draft of the new CJCSM 3150.23 Joint 
Reporting Structure (LOGFACREP) to improve the process of Service inputs to 
the LFF. 

@ Establishing a routine process for maintaining the LFF. 

e@ Updating supporting files, such as the TUCHA. 


@ Developing a procedure to use DoD class/subclass supply codes. 


The intent of this thesis is to devise a method of analyzing historical data for Navy 
surface ships by fitting a mixture of probability distributions to this data. The use of these 
distributions should provide accurate consumption and sustainment rates for these ships. 
Included in these results will be the expected values and variances of the respective 


sustainment rates. From these distributions, notional ARGs and CVBGs can be assembled 
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for planning purposes, and the sustainment rates aggregated for the entire BG. 
This thesis will be narrowed in scope to model only the consumption rates for the 
seven subclasses of Class LX, Repair Parts. It will also be limited to four classes of ships that 


might compromise a typical CVBG. These ship types include: 


@ CVN-68 (Nimitz Class) Aircraft Carrier, 


® CG-47 (Ticonderoga Class) Guided Missile Cruiser, 


DD-963 (Spruance Class) Destroyer, 


FFG-7 (Oliver Hazard Perry Class) Guided Missile Frigate. 


Newer classes of ships (such as DDG-51 Guided Missile Destroyers) were not included due 
to lack of data, and several older classes of ships (such as FF-1052, CG-16 and DDG-2) were 


excluded, despite a wealth of existing data, because they are being decommissioned. 
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IV. DATA COLLECTION AND ANALYSIS 


The data used to model the consumption rates of repair parts for this thesis was 
provided by SPCC Mechanicsburg, PA. SPCC is one of the two Inventory Control Points 
(ICPs) within the Navy Supply System, and it is responsible for managing all non-aviation 
related items in the Navy’s inventory. The data consisted of all surface ship requisitions from 
1 January 1992 (Julian date 2001) to 29 September 1994 (Julian date 4272). There were a 
| total of 7,020,303 requisitions. 

Additional information provided by SPCC was weight and cube data on items in both 
the Navy and DLA supply systems. The records in these fields are listed by NIIN. There are 
469,361 records in the Navy file and 3,533,161 records in the DLA file. 

A. DATA MANIPULATION 

The data files provided were stored on 6250 bpi (bits per inch) magnetic tapes in fixed 
block format. For this reason, plus the sheer magnitude of data, the files were loaded onto 
the Naval Postgraduate School (NPS) mainframe computer. All of the data manipulation that 
follows in this section was performed with SAS code or in VM/CMS with the assistance of 
consultants at the NPS computer center. 

The first step in manipulating the data into usable format was eliminating requisitions 
from all ship types except the four classes of ships previously mentioned. This process 
actually involved several steps, the first of which was determining the desired UTCs from the 
TUCHA. The second step was to find the Unit Identification Codes (UICs) of the individual 


ships within the desired UTCs. This was accomplished using a UTC to UIC conversion Table 
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also provided by SPCC. Appendix A lists the four desired UTCs and their respective UICs. 
Code was then written in SAS to sort through all of the requisitions, eliminating those 
with undesired UICs. From the requisitions that remained, only the following fields were 
saved (the number in parentheses represents the number of characters in that field): NIIN (9), 
UIC (5), Julian date of requisition (4), quantity ordered (5), FSG (2), COG (2), and a 
deployment indicator (1). The FSG and COG are supply codes that provide additional 
information about the items in the supply system, such as the manager of that item. The 
reason for copying these fields will be discussed later. The po ee indicator is a binary 
character, with a one being used if the unit submitting the watiies was deployed, zero 
otherwise. As a result of eliminating requisitions from ship types that we are not modeling, 
the following numbers of requisitions remained: 358,918 (CVN-68); 276,386 (FFG-7); 
247,816 (DD-963); and 213,615 (CG-47). | 
The next major step involved assigning a weight to each of the requisitions. However, 
after combining the DLA and Navy weight files and eliminating duplicate entries, 
approximately 40% of the 3,606,482 unique records did not have weights assigned. It was 
necessary to devise a policy for adjusting to this lacuna. Merely calculating the average of 
all the records and using this value in place of all the missing weights was ruled out because 
this would significantly alter the resulting distributions and lead to underestimation of the 
variance. Rather, the records were broken into two categories, those with weights and those 
without. Random sampling (without replacement) was performed from the list of NIINs that 
had recorded weights, and these random values were inserted into a field called weight. This 


method of sampling generated a random permutation of weights which was used serially 
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whenever a weight was needed. In this way, the same part could have different weights 
assigned in subsequent draws. Otherwise, the correct weight of the item being requisitioned 
was inserted into this field. 

In an effort to reduce the amount of random sampling required, as well as improving 
the accuracy of the sampling, we limited our random draws to only NIINs that had been 
requisitioned during the three year period. This approach prevents the chance that an 
extremely heavy item carried by DLA for another service could have its weight inserted into 
a Navy requisition. This technique, plus the previously discussed method of sampling each 
time we came across a NIIN without a weight, was deemed the best way to minimize 
distortion of the original data. Out of the 1,096,735 requisitions, there were only 174,225 
unique NIINs requisitioned. From these unique NIINs, there were 135,970 (78%) that had 
weights assigned, while only 38,255 (22%) had missing values. Each occasion where a 
requisition occurred with a NIIN that had a weight of zero, a random sampling from the list 
with weights re a value into the weight field of that requisition. 


Table 3 shows the results of the above procedure for each of the four ship types: 


. requisitions with weight weight w/out weight | w/out weight 
cea | risers] sass] na] agooa | 89 
; : 6.2 















an at ae a ad 
358,918 | 304,458 $4,460 


Table 3. Results of Random Assignment of Weights to NIINs with Missing Values. 


The results achieved by using this methodology are significantly better than using the average 
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values of all 3.5 million items (40% of which do not have values). For each requisition, the 
fields of quantity and weight were multiplied together to achieve a total requisition weight. 

The third major step also involved a very subjective methodology that is worth 
describing in detail because it too forms an underlying assumption for the entire model. This 
step sorted through each of the requisitions (for each of the UTCs) and discarded those that 
did not belong to Class [X (Repair Parts). Those requisitions identified as being from Class 
IX were broken down into one of the seven subclasses listed in Table 2. These procedures 
utilized a table which converted the FSG and COG into a corresponding DoD supply class. 
This table is included as Table 4. This Table is subjective in that every item in a FSG/COG 
grouping is assigned to the same DoD supply class. This is not what the FSG/COG codes 
were designed for, so individual items may wind up being assigned to a subclass to which they 
do not belong. However, this method does allow generalized groupings of supply items into 


class/subclasses, allowing the Navy to comply with Joint Staff directives. 
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FSG COG DoD Supply Class 


a ee ee) eS 


Table 4. Guidelines for Assigning Repair Parts to DoD Supply Classes. 





The data was now broken down into 28 separate files (four UTCs, seven subclasses 
of repair parts per UTC). Within each of these files, the weights from multiple requisitions 
by the same unit on the same date were added, resulting in a single value. The final step in 
preparing the data for analysis was to use a spreadsheet to edit out days of requisitions from 
units that were not deployed. This was also a fairly subjective step, because different ships 
reported differently as to whether they were deployed or not. For instance, several ships 
considered five day underway periods as deployed, and a few ships reported themselves as 
deployed for 1000 consecutive days. Other interpretations of whether or not a unit was 
deployed were involved in this step, but consistency was used throughout. After this was 
accomplished, the probability that an item from each subclass of Class IX was calculated for 
each of the four UTCs. These results, along with breakdowns of the percentages of 
requisitions involving actual weights versus sampled weights is included as Appendix C. 

B. DATA ANALYSIS 

The original intent of this thesis was to fit a single probability distribution to the 
weights of the daily requisitions in each subclass of repair parts. A mixed distribution would 
result, with a probability p, of having no requisitions on a given day (hence zero pounds being 
requisitioned), and another probability p, of a eauianon occurring on that day. The weight 


of that requisition would be described by the probability distribution that best fit that set of 
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data. The resulting expected daily requisition weight (for a given subclass of repair part and 


ship type) would be computed: 


E[weight] = ip, x 0} + p, x fine) [1] 
0 


where f(x) represents the probability density of the daily positive requisition weight. In 
addition to the calculating this expected value, we could also calculate its variance. 

From the beginning of the study, it became apparent that this mopisucl was not 
appropriate for Class [X(D), administrative vehicles. Regardless of ship type, the number of 
requisitions for this type of repair part are so low (a maximum of 12 requisitions during a 
three year period) that 1t was decided to model this subclass using weighted averages. 
Equation 21 in Chapter V illustrates the modifications to Equation 1 to arrive at an estimate 
for this subclass. However, Equation 1 could be not be applied directly to most of the otha: 
data files either (exceptions being Class [X(L) for FFG-7 and DD-963 ship types). The 
reason this simple node will not work for the other data sets is due to difficulties in fitting. 
The Weibull values are polarized because it is impossible to adequately fit a distribution to 
them. This is due to the fact that most of the daily requisition weights in any subclass for any 
of the four ship types were very light (ten pounds or less). On the other hand, there = 
extreme cases of requisitions whose weight exceeds 100,000 pounds in almost every data set. 
Moreover, these outliers are from real requisitions which, far more often than not, are from 
actual weights and not the randomly sampled weights substituted for the NIINs listing 
weights of zero. 


Figures 3 and 4 on the following page illustrate how these extreme weights prevent 
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Figure 4. Histogram of CVN-68 Class IX(A) Data and Weibull(.47, 114) 
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any good fit of a distribution to the data. The data in the figures is from CVN-68 (aircraft 
carriers), Class [X(A), and it is representative of most of the data sets. Figure 3 is a plot of 
the empirical cumulative distribution function (ECDF) of the daily requisition weights, while 
Figure 4 shows a histogram of the same data and an attempt to fit a Weibull distribution to 
that histogram. The extreme weights force most of the data into the first class, or ‘bin’, of 
the histogram, preventing a good fit. Attempting ad hoc adjustments with the number of 
classes in the histogram did nothing to relieve this problem. 

1. Extreme-Value Distribution 

The method chosen to handle the cases of the outliers is that of modeling them using 
one of the extreme-value distributions for maximums. Since weights are positive, the 
extreme-value distribution function chosen for our data sets is the Type II, or Cauchy type, 
which is defined as: | 

G(y) = exp(-y") for y>0, y>0. [2] 

This form is derived as the limiting distribution of a normalized set of maxima 
Y(n) = a [Ref 14]. This type is best suited for maximums from densities having 
thick tails, aa as the extreme values in the data sets of the requisition weights. 

If these very heavy (and rare) requisitions can be modeled using an appropriate 


distribution, the resulting daily expected requisition weight becomes the mixture: 
E|weight] = (p, x 0) + (By X Hy) + 3 * Hy) [3] 


where p, is the expected value of the “outlier free” portion of data, and p, is the mean of 
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those observations in the extreme.’ The values p,, p, and p, are the mixing probabilities, and 
are defined as: p, for the days without requisitions, p, for the requisition weights not 
containing any outliers, and p; for requisitions of extreme weights. 

Figures 5 and 6 on the following two pages illustrate the effects of removing the 
extreme values from the data set representing CVN-68 ships, Class IX(A). After truncation 
of the 19 heaviest requisitions (which will be modeled with the extreme-value distribution), 
the remaining data yields a smoother ECDF and a histogram which appears better suited for 
being fit by a common probability distribution. 

2. Weibull Distribution 

A method was devised to identify the extreme values for all of the data sets except 
those for Class [X(D) (all ship types) and Class [X(L) (FFG-7 and DD-963 ship types only). 
The extreme points were removed following this procedure, and after this truncation, the 
resulting ECDFs appear to resemble that of a Weibull distribution. Weibull distributions are 
‘J-shaped, ai provide a richer set of models than the Exponential type. They are often used 
in probability when attempting to model the time required to complete a task or the time to 
failure of a piece of equipment [Ref 15]. The data in the tail of the remaining requisitions was 
another factor in the selection of the Weibull distribution. 

The Weibull distribution describes a continuous random variable X with two 


parameters: 0 and B. The range for both of these parameters is 8,8 > 0. The probability 


' The use of the Type II extreme-value distribution does not contribute a finite p, 
for Equation 3, but we shall devise a method of dealing with this issue later. 
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Figure 6. Histograms of data sets for CVN-68, Class IX(A), as extreme values 
are removed. Y-axis represents number of observations. Note changes in the 
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density function (pdf) for a Weibull distribution 1s of the form: 


x\b 
B) me if x>0 (4) 


F(x;9,B) = (4 
0 otherwise 
The parameter B is called the shape parameter, while 0 is the scale parameter. The 


Cumulative Distribution Function (CDF) of a Weibull distribution 1s: 


x\P 
F(x;0,8) = 1-e ti) for x>0 [5] 

The parameters 6 and 8 for the Weibull distributions will be estimated using the 
method of maximum likelihood. The likelihood function can be described as “giving the 
probability of observing the data as a function of the parameters. The Maximum Likelihood 
Estimate (MLE) of a set of parameters is that value of the set that maximizes the likelihood - 
that is, makes the observed data ‘most probable’ or ‘most likely’ “ [Ref 16]. The derivation 
of the MLEs for the Weibull follows over the next few pages. 

Recall the density function for a Weibull distribution from Equation 4. Let the n 
values from any given data set be x,, x,, ..., x,. Assuming that x,, ... , x, are independent and 
identically distributed according to a Weibull(6, B) distribution, then their joint density is the 


product of their marginal densities: 


x\P 
f(,,%,...x,|9,B) = [] ES xP le | [6] 


i=] 


The likelihood function (L), when regarded as a function of 6 and 8, is: 
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n [_7 E yy x; \p 
(0,8) = IT x] (a) 


i=] 


The log likelihood function LL = In[ L(6,6) J, is: 


n n (x.\B 
LL = nin(B) - nBln®) + (B-I)y Ine) - (*) 


f=] 


The partial derivatives of LL, with respect to the parameters 8 and 8, are: 
OLL -nB . ( *) po sf 
ee ee ws = | eee 
36 0 s P 60 on Q2 


OLL 


8 pM) + > InGs) - 3 & | nf = 


Equations 9a and 9b are set equal to zero, and after some algebra, yield: 


—J 





Dees 
eh = Ll 
n 
d xPin@) SY ine) 
- J=1 _ dsl 
n n 
Dr, 
i=] 
The first parameter we solve for is B. We do this by defining g(B), where: 
> xP inc) 
g(B) = ——_—_ - — 
d xP 
i=] 
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An iterative process (Newton’s Method) is now used to solve for f: 


y In(x,) 
B, = B, - e(8) - ——— [12] 
g(B) 


where g‘(B) is the derivative of g(B): 


n 2 
> xP nt) 
i=1 








[y. :!) bs xP In(x,)? 


fh) ce iii 
baal 


B, is initially set equal to 1, and replaced by 61 at each iteration, and the process of updating 
B, continues until |B, - B,| < €, where € is a very small number (we use 0.001). Once the 
value of B has been determined using this method, the value of 0 can be calculated using 


Equation 10a. 


Having calculated these estimates, their values will be inserted in the following 


equations to get summary statistics. The mean (E[X]) and variance (Var[X]) can be 


determined from the formulas: 
E(X) = ar 1 + x) 14] 


Var(X) = ert ' 2) : rf + 4) [15] 
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where I(x) denotes the gamma function, which is defined for all real x>0 by: 
T(x) = [ t* leat [16] 
0 
Finally, the 100 x p™ percentile of a Weibull distribution is of the form: 
i 
X, = @|-In(1 - p)/? [17a] 


We will also be working with the median value of these Weibull distributions, which is defined 


as the 50" percentile (or X, = .5). This special form of Equation 17a is: 
- | 
te @|-In(.5)]? [17b} 


3. | Methodology for Determining Extreme Values 

We wanted to develop a standard methodology for identifying extreme weights which 
could be used on all of the data sets, rather than using subjective guesswork. Such a universal 
standard has eluded us, nor was it possible to identify a single weight as the cutoff between 
regular and extreme. However, we were able to establish rough “rule of thumb” guidelines 
that enabled us to set criteria which could generally be applied to a given subclass of repair 
parts (across all ship types). 

The first step in this process involved calculating the ratio of the empirical mean and 


median values from a given data set. These are compared with the same ratio using the 
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theoretical mean and median values. Such ratios provide information about the skewness of 
the underlying data sets. Upon completion of this step, the largest (1.e., heaviest) item in the 
data set 1s removed, and the calculations repeated until we get a reasonable match between 
the two. The idea here is to remove outliers until the empirical ratio of mean to median 
agrees with the theoretical. This process generally repeats around 50 times; the actual counts 
ranging from 11 to 85. Operationally, the ratios of actual mean to actual median and 
theoretical mean to theoretical median are simultaneously plotted against the percentage of 
Tequisitions being truncated. Figure 7 on the following page illustrates the results of these 
plots for Class [X(G) repair parts for all four classes of ship. 

Our original intentions were to examine these plots and identify obvious “knees” in 
the curves of the empirical skewness (1.e., when the ratio of mean to median stopped dropping 
off sharply as additional weights are truncated). These knees in the curve would serve as 
appropriate places to stop truncating requisition weights as extreme. This method only 
worked in a few cases, so we retreated to the method of the previous paragraph. The 
intersection points are not sharply determined, but approximate solutions seem serviceable. 
This approach worked well for some of the subclasses (such as IX(G) in Figure 7), but not 
in all of them. In several cases the plots of the ratios behaved very unpredictably. 

We then sought ways of developing criteria allowing for the truncation of extreme 
points from the other data sets. In doing so, we developed the following rough guidelines: 

@ As the actual ratio of mean to median decreases to approximately 2.5, the behavior 

of the skewness curves starts to behave unpredictably. The reason for this is 


mainly due to the fact that the values being truncated are no longer extreme. This 
type of activity is depicted in Figure 8 on page 38: 
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Figure 7. Plots of actual and theoretical ratios of mean to median against 
percentage of truncated requisitions. Class IX(G) is plotted for all ship types. 
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® For the larger data sets (Classes IX(B), [X(G), [X(T)), do not truncate more than 
5% of the weights as extreme. The remaining data sets (Classes IX(A), [X(L), 
IX(M)) should not be truncated beyond 13% of the weights. 


FFG-7, Class IX(M) 


Percentage of Truncated Requisitions 





Figure 8. Illustration of Low Ratios Between Mean to Median 
Values. 


@ Calculate the ratio of the actual (mean to median) ratio over the theoretical (mean 
to median) ratio. Continue truncating weights as long as this value is greater than 
2.2. Further truncation will probably be required, but this serves as a good upper 
bound. 


@ Lastly, observe the values of the Kolmogorov-Smirnov statistics as additional 
weights are truncated to see if a better fit of a distribution to the data can be 
achieved by further truncations. A discussion of this statistic follows in the next 
section. 


These guidelines should not serve as rigid rules to be followed for all circumstances. Rather, 


they serve as steps to point future analysis in the right direction. 
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4. Kolmogorov-Smirnov Statistic 
The Kolmogorov-Smirnov (K-S) statistic measures the goodness of fit of a 
hypothesized CDF to a given sample ECDF. In other words, the K-S test evaluates how well 
our Weibull(0,8) distributions fit the ECDFs of the data sets. It accomplishes this by 
determining the maximum difference between the sample ECDF (the requisition data) and the 
hypothesized CDF (the Weibull distribution) [Ref 14]. The K-S distance function (D,) is 


defined for a given H,: X~F(X) as follows: 


+ Max l 2 | 
Dn = Veien ‘4 Xo} nee] 
- _ max _-t) 
Dn > Veicn HX) aD eu) 
D,= max|D Dt [18c] 


where X,1, X(, ... , Xqy are the n order statistics from the data set, i is an index from | to n, 
and F(X) is our hypothesized distribution. 

The result of the K-S test is the distance function D,, and small values indicate better 
fits. “The test rejects the null hypothesis (H,) if D, exceeds some constant d, ,.,, where a is 
the specified level of the test” [Ref 15]. Our calculations used a value of a=.05. The K-S test 
was selected over other more common tests (such as the chi-square test) because it is easier 
to use on large sample sizes and it uses ungrouped data, so the goodness of fit is measured 


at every data point. It also prevents the necessity of having to determine interval widths and 
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does not require a correction for continuity [Ref 17]. 
The K-S distance function is evaluated only on the data with positive (i.e., greater than 
zero) values. Recall from Equation 3: 
P, = probability of no requisition on a given day 
Pp, = probability of a requisition of non-extreme weight 


Pp, = probability of a requisition of extreme weight. 





We next define: 
n, = number of requisitions described by Weibull 
n, = number of requisitions described by EV function 
OST De cece Vt 
P 
Ps eee ae 
1 -p, 
P 
Pp : = ; ? 
I - p, 


and let Xq, Xo, .-., X2+n3) be the ordered weights of the requisitions in any given data set, 
where Xq) < Xi < ... $ Xqo4n3, We further define F(X) as the Weibull distribution and G(X) 
as the EV distribution. Now we calculate F,(Xj) for all i= 1, 2, ... , m+n, where F,(X,)) is 


defined as: 


F(X) = (PixFX)} + psx GX) [19] 


The value of F,(X,)) is used in Equation 18c to obtain the K-S distance of our function. 
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De Fitting Distributions to the Data 


Table 5 lists the breakdown of the data into the three categories of requisitions: 


Class | Days No | Days Days 

Deplo Req | ~Weib ~EV 
Ix(A) eee 
cw | 1x®) ee eee he! 


| os] 02 

xc) | 918 | 801 00 

IX(M) | 918} 708 001 
sts | 356 


IX(T) 


067 | 
a 
am 
387 
rx(a) | _ 6364] 6216 | 977 | 
coa7 | 1x@e)| 6364 | 3831 | 602 
1x |_ 6364] 2957] 465, 
xq) | 6364] 6319] 992, 
rx) | 6364] 6098 | 957 
caa7 | ix | 6364] s209| 819. 
FrG7 | 1x(a)| 407] sat | 97 
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548, 
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Class Days a | Days i oe 

Deplo | ~Weib ~EV 
z= xem) |_s344| 7986 | ns} a8} ots 
g44| 7136] 855] 7 | tat 


Lo. 5. Results of Breaking Data into Three Gatesories of at] 





No distributions are being fit to the data from Class IX(D) due to the very low 
numbers involved. From the remaining data sets, most have been broken down into two 
subsets: the ‘regular’ requisition weights that will be fit by a Weibull distribution, and the 
extreme cases, which will be handled by the EV distribution. A small number of remaining 
data sets will be modeled only with the Weibull distribution (i.e., p, is zero). 

Determining the parameters for the Weibull distribution is straightforward, using the 
method of maximum likelihood previously discussed. The parameters of the Weibull 
distributions fit to the data sets are listed in Appendix C. Included in this Appendix are the 
old K-S distances determined prior to the truncation of extreme points, as well as the K-S 
distances calculated after the truncation. Except for the two instances where no truncations 
of extreme points occurred, the fit of a Weibull distribution to the data is better. In some 
cases, this improvement is quite dramatic. There is difficulty in determining the parameters 
for the EV distribution. Recall from Equation 2 that G(y) is the limiting distribution of a 
standardized Y = ae where X is the requisition weight and a,b are scaling parameters. 
The two parameters a and b are both unknown, as is the value for y in Equation 2. These 


values will be estimated using an ad hoc technique. 


The technique applied is that of adjusting the parameters a, b and y until the Q-Q plot 
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(empirical versus theoretical) of the data is as close to a straight line as possible. We 


accomplish this by plotting v; against u; for each data set, where v; and u;, are: 





v= nf -n(—£ ) [20a] 
u, = In{X, - B), [20b] 


where m 1s the number of values in the data set; j is a vector of indices from 1 to m. Next 
we adjust the value of b in Equation 20b until the scatterplot foams as straight a line as 
possible. Simple regression is performed on this line to obtain the slope and intercept. A 
demonstration of this technique is provided in Appendix D. Submitted without explanation, 
the slope of this line is the negative value of y. The value of a is calculated from the 
intercept, which is equal to yln(a). The results of this technique are also listed in Appendix 


D. 
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V. RESULTS OF ANALYSIS 


A. DATA SETS TREATED AS WEIGHTED AVERAGES 
The first section of this chapter deals with the very small data sets of Class IX(D). 
These data sets had very few requisitions, and those that did occur had relatively low 
weights. We will call the number in the last column the Central Value. This Central Value 
is calculated with only a slight modification to Equation 1: 
E| Central Value | = (p,x0) + (—D,*p") | [21] 
where 1” is the empirical average requisition weight from all of the requisitions in that data 


set. The values for this subclass are extremely low, as Table 6 illustrates: 


Deployed | With Req’s | Given Da Weight Value 


Table 6. Expected Daily Requisition Weight for Class [X(D). 















B. DATA SETS WITH NO EXTREME VALUES 

This section discusses results of the cases where no requisition weights were 
truncated. These were both in Class IX(L), for FFG-7 and DD-963 ship types. The 
parameters of the Weibull distribution fitting those data sets are given in Table 7, as well as 


the resulting Central Value, which is solved using Equation 1 and the mean value of the 
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Weibull distribution that was fitted to the actual requisition weights. The K-S distance of 
the distribution fit and the rejection cutoff for the K-S test at a value of a=.05 are also 


included: 


Sl 
: Value dist cutoff 
Ce ee eee ee 
pps63 | _.o10| 1623] _.73| 19.772 | 760.08 | 198] _.130| _.152. 


Table 7. Expected Daily Requisition Weights for Data Sets with No Truncation. 






C. DATA SETS WITH EXTREME VALUES 

This section discusses the results of most of the data sets that were modeled, and it 
includes data sets which include days with no requisitions, days with requisitions of items of 
‘regular’ weight that are modeled by the Weibull distribution, and days with extreme 
requisition weights which are modeled by the Extreme-Value distribution. The Central 
Values for this section were calculated using a modification of Equation 3. The point that 
must be addressed is that no mean and variance exist for the Type II Limiting Distribution 
that we chose to use. Instead, we calculated the three quartiles of the distribution. That is, 
Equation 22 is used to calculate quantiles for values of p = .25, .50 and .75: 

p =e -y) | [22] 

The quantile at p=.5 (or q,5) ) is the median value that we will use, after scale change, 
in place of for the purposes of calculating the Central Value. No variance exists for this 
distribution either, and we will calculate the Inner-Quartile Range (IQR) for use as a measure 


of spread. It is defined as: 
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TOR = 475) ~ 425) [23] 
The IQR gives an interval of values that contains 50% of the cases. The value y in Equation 
22 comes from Table D-1 in Appendix D, and the equation is solved for y. The solution to 
Equation 22 must be converted to the original units for use in the Tables using Equation 24: 

X=b+aY [24] 
where the values for ‘a’ and ‘b’ come from Table D-1. Thus, q,; is converted to the 
“xX” scale and used for pn, 

_ The effect of all this is to output some hybrid measure. The choice of a central value 
is awkward, and requires further study. The expected value is preferable because the period 
total is computed as the product of expected value with days in the period. The median (or 
other quantile figures) does not have this very desirable property. Indeed, for many of our 
cases, p, > 0.5, and for these instances, the median is zero. As an expedient, we are using 
the Weibull mean for 1, and the extreme value median for p,. At least this latter value is 
finite and serves as something that can be used in computation, albeit we anticipate the result 
to be on the low side. With regards to variablity, one can be reminded that for normal 
variates the IQR is roughly 1.35 standard deviations. This conversion 1s of limited usefulness 
for heavily skewed distributions. Although the entries do not blend well, they are the best 
we can do under the circumstances. 

The output in Table 8 on the next page illustrates the results of applying the 
preceding method. 
A brief discussion of what the output in Table 8 represents follows: 


Column 1 - The ship type and Subclass of Repair Part being modeled; 
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Columns 2,3,4 - The probabilities defined in Equation 3; 

Column 5 - Expected Values of the Weibull distributions fitted to the data; 
Column 6 - Product of p, and the values from Column 5; 

Column 7 - Variance of values from the Weibull distributions; 

Column 8 - The q,;) quantile (median value) from E-V distribution; 
Column 9 - Product of p; and the values from Column 8; 


Column 10 - The IQR of the extreme weights (q.75) - 4¢25)); 






Class [pt {p2__[p3__|u,__fp2*u,__|VarfWeib] |q(.5) __{p3*q(.5) IQR 
CV 9A {0.36 |0.64 |0.00 |222 [142.08 [246420 _ [100352 __|100.35 


Cv 9G_|0.07 [0.86 [0.07 [411.82 [355.81 [722656 __|113856.5 |7856.1 












cv 9m [0.77 {0.23 0.00 [83.08 |18.94 [25423 |soo00__—go. [N/A 
cv 9T 0.39 |0.58 |0.03 [227.33 |132.08 [375967 _ |28417.8 
icG9A_|o.98 [0.02 [0.00 [3.76 |0.08 [15.64 (87.35 —|o.17__—fo0.89 
cG 9B |0.60 |o.40 |0.00 {56.16 [22.24 |10616.6 _|19734.6 _|39.47 
cG9L_|0.99 [0.01 |0.00 |12.71 |0.09 [179.16 _—fa50__——0.25 
cG9M 0.96 |0.04 |o.00 |40.51 [1.7 [3397.2 390__—ifo.39 
cG9T_|o.82 [0.17 |0.01 [26.21 [4.56 [1883.7 [2149.9 [15.05 
FFG 9A [0.98 |0.02 [0.00 [4.7 [0.1 (33.28 _—|611.18 
FFG 9B |0.67 |0.33 {0.00 [54.57 [17.74 [12687___|15786.4_|31.57 
FFG 9G |0.55 {0.45 [0.01 [47.15 [20.98 
FFG 9M|0.96 |0.03 [0.01 [12.11 [0.4 [162.5 (118.79 |0.59_—[63.16 
FFG 9T |0.86 [0.14 |0.01 [22.6 
DD 9A_|o.98 |0.02 |o.00 [4.52 0.07 [24.72 __—*{75.68 _—*|0.08_—(59.49 
pp 9B {0.65 |0.35 |0.00 155.61 |19.57 |10409.25 |7495.3  |14.99 _ |6910.9 

DD 9M [0.96 |0.04 j0.01 [14.02 10.53 {249.12 [132.21 
Dp 9T_|o.86 |0.14 |o.00 [26.95 [3.8 [2444.1 foss.8 3.84 [2215.04 _| 





Table 8. Parameters for Use With Mixed Distributions. 
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D. COMPARISON OF OLD LFF CONSUMPTION RATES AND NEW RATES 
Table 9 below illustrates the significant differences that exist between the old factors 
and the new factors calculated using the mixed distribution methodology developed in this 


study: 


Current Current Current 
CVN-68 CG-47 


Factor Factor Factor Factor Factor Factor 
78.00 616.36 61.71 34.57 


0.14 
0.02 
0.15 


463.95 





Table 9. aera of Carre Factors j in Navy LFF To New Sache Derived 
from Recent Historical Data and Mixed Probability Distributions. 

This table demonstrates that the current factors are generally low by a factor of ten. 
This could mean that existing OPLANs are deficient in estimating sustainment requirements 
for a CVBG by the same factor of ten, which calls into question any logistics feasibility 


analysis that has been performed on those OPLANSs. 
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VI. CONCLUSIONS 


It is the conclusion of this thesis that the Navy’s LFF is grossly out of date. Using 
available historical data, we have justified a rough rule that the current consumption factors 
in the LFF are “low by an order of magnitude”, at least in the case of Class IX, Repair Parts. 
This conclusion is underscored by the fact that the data was collected during “peacetime” 
deployments. The operational tempo of these deployments may have been high, but the ships 
were not attacked by the enemy. Logic dictates that once battle damage occurs, the requests 
for repair parts will only increase. 

A significant amount of data massaging was required by this study. Several ad hoc 
techniques were necessary because of the condition of the data. While we submit that these 
breaches of pure statistical analysis wee carried through the analysis - the results, at jest 
in form, they cannot explain the order of magnitude difference between the old and new sets 
of factors. However, unlike the current factors, we developed and clearly documented our 
methodology. As a result, future analysts in this area can use all, part or none of our 
methodology, because they will know what we did to arrive at our results. 

The Navy is moving to comply with all aspects of the directives issued by the Joint 
Staff regarding the improvement and maintenance of the LFF. One of the most significant 
results of the work done by PRC, Inc., was to identify specific organizations to be responsible 
for maintaining the factors current for individual subclasses of supply. Additionally, a single 
person has been established as a point of contact at each of the cognizant agencies. The 


results of these steps should make for easier and more efficient maintenance of the LFF. 
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The requisition data that was provided by SPCC for this study was complete and 
thorough. We feel strongly that historical data such as this 1s the best way to model and 
predict future consumption factors. However, weight and cube (volume) data must be 
accurate to successfully model the consumption factor as a function of weight. We 
recommend that any further studies in this area begin by cleaning up this data. Once this has 
been done, it can be followed by a more studied and carefully structured approach to the 


model fitting. 
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APPENDIX A. UNIT TYPE CODES AND UNIT IDENTIFICATION CODES 


This appendix lists the UTCs of the four ship classes being modeled, as well as the 


UICs of the ships belonging to that UTC: 


UTC 5CVN1 (8 CVN-68 Nimitz Class Aircraft Carriers 








03368 
03369 


20993 
21247 


21297 
21412 


21847 
21853 


UTC S5DDD3__(31 DD-963 Spruance Class Destroyers) 


20574 20591 20611 20834 20589 
20575 20598 20612 20835 20590 
20576 20599 20613 20836 20603 
20586 20600 20614 20837 30604 
20587 20601 20615 20838 20617 
20588 20602 20616 20839 20833 
21416 


UTC 5GEGI1 (52 FFG-7 Oliver Hazard Perry Class Frigates) 


21028 20975 21103 21231 20973 
21032 20976 21104 21232 20974 
21033 20977 21105 21233 21058 
21034 20978 21106 21234 21059 
20964 20979 21107 21235 21200 
20965 21052 21108 21236 21201 
20966 21053 21109 21350 21430 
20967 21054 21110 2135] 21111 
20968 21055 21197 21352 : 
20969 21056 21198 21390 

20972 21057 21199 21391 


UTC 5JCG2 (27 Ticonderoga Class Aegis Cruisers) 


21281 21387 21450 21658 21346 
21225 21388 21451 21684 21449 
21295 21389 21623 21827 21657 
21296 21428 21624 21828 
21344 21429 21625 21829 
21345 21656 21830 


21447 
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APPENDIX B. DATA MANIPULATION RESULTS 


The four tables in Appendix B display the results of the data manipulation discussed 
in Chapter IV. They show the number of days deployed units within each of the UTCs 
requisitioned a repair part from each subclass. Additional information in the tables include 
the number of requisitions that had actual weights assigned and the total number of 
requisitions. From this information, the percentages of requisitions based on actual weight 


data is compared to the percentages of those that had random weights assigned. 


Total Percent Percent 

of reqs w/ | number reqs =| reqs w/ 
actual of reqs | w/actual | random 
weight weights weights 


t 
7 





Table B-1. Result of Data Manipulation on CVN-68 Class. 
918 Deployed Days. 
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Number | Number | Number | Percent Percent 

of days | of reqs w/ | of reqs w/ reqs reqs w/ 

with reqs | actual random | w/actual | random 
weight weight 






179 

cx | 3) aps 
rae) | aor] somo] aimes| ois] 82 
rxay | 8] ml 1o00) 00 


Table B-2. Results of Data Manipulation on CG-47 Class. 
6364 Deployed Days 







| ts 
1x(A) 199} 900] 100 
3 


Number | Number | Number | Percent Percent 

of days | of reqs w/ | of reqs w/ reqs reqs w/ 

with reqs actual random | w/actual | random 
weig 





Table B-3. Results of Data Manipulation on FFG-7 Class. 
8,407 Deployed Days. 
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Number | Number | Number | Percent | Percent 

of days | of reqs w/ | of reqs w/ reqs reqs w/ 

with reqs actual random | w/actual | random 
weight weight weights 







rma |e) i) m) oar 
xq) | 2081] aire) aaars| our 
rixe) | sof sof a] ore 
xn | 1208] 3579] 37m | _ 949 


Table B-4. Results of Data Manipulation on DD-963 Class. 
8,344 Deployed Days. 


dl 


weights 
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APPENDIX C. RESULTS OF FITTING WEIBULL DISTRIBUTIONS. 


Appendix C lists the parameters of the Weibull distributions used to fit the data sets 
after the truncation of extreme values. It also illustrates the improvements in this fit by listing 


the old K-S distance prior to truncation as well as the new K-S distance: 


New K-S 
Distance 


No Trunc 


° 


076 
148 
062 
068 
137 
080 
076 
115 
079 












Ra Ma ON a 9 
| Distance Distance _ 
Tops | mw) | se | asi] ow | ose 
Pops | meg |e | sa) atom 


Table C-1. Results of Fitting Weibull Distributions To Data Sets 
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APPENDIX D. DETERMINING PARAMETERS FOR EV DISTRIBUTION 


This appendix provides an example of the ad hoc technique used to determine the 
parameters for the EV Distribution. Figure D-1 below illustrates the effects of varying the 
parameter b in equation 19b to form a straight line in the scatterplot. The example shown is 


for FFG-7 ship type, Class IX(T). The value of b=200 yields the straightest line: 


FFG-7, Class IX(T) 





Figure D-1. Effects of varying parameter b to get a straight line. 


Once the best value of b has been determined, the slope and intercept of the line are 
calculated. These values are used to find the values of the other two parameters in the EV 


distribution. A linear regression performed on the data provides the desired values, as is 
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demonstrated in Figure D-2: 


FFG-7, Class IX(T), b=200 


=J 





Figure D-2. Calculating slope and intercept for EV distribution. 


Table D-1 displays all of the parameters for the EV distribution that were obtained by 


using this technique: 







[eves | sol ss[ save us| as 


CV IX(T) 3000] —-0.7368 0.7368 | 15,456.20 
CG IX(A) Pg | -1,5918 6.6741 1.5918 66.208 
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| Class |b | Slope __| Intercept | yy | a 


CG IX(B) ~ 4000 -1.1409 | 10.6617 1.1409 | 11,411.37 
CG IX(G) 3000 -1.0455 10.6797 1.0455 | 27,307.62 












a Oe 
= a 
— 
ae 
oe eT 
oe eC a ae 
= 
DD IX(M) 

DD IX(T 0.6269] 3.58831 0.6209 





Table D-1. Parameters for EV Distributions Obtained by Ad Hoc Technique. 
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